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It is shown that superpositions of path integrals with arbitrary Hamilto- 
nians and different scaling parameters v ("variances") obey the Chapman- 
Kolmogorov relation for Markovian processes if and only if the corresponding 
smearing distributions for v have a specific functional form. Ensuing "smear- 
ing" distributions substantially simplify the coupled system of Fokker-Planck 
equations for smeared and un-smeared conditional probabilities. Simple appli- 
cation in financial models with stochastic volatility is presented. 
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1. Introduction 

One often encounters in practical applications probabilities formulated as 
a superposition (or "smearing" ) of path integrals (PI) of the form 

P(x b , t b \ Xa , t a ) = / dv LJ(V, t ba ) / VxVp dr(«p4-H(p,x)) _ (1) 

JO Jx(t a )=x a 

Here tu(v,t ba ) is some positive, continuous and normalizable function on 
K + x IR + with t ba = tb — t a being the time difference. Examples of (1) can 
be found in financial markets, 1,2 in polymer physics, 1,3 in superstatistics, 4,5 
etc. 

Whenever a smeared PI fulfills the Chapman-Kolmogorov equation 
(CKE) for continuous Markovian processes, the Feynman-Kac formula 
guarantees that such a superposition itself can be written as PI, i.e. 

P(x b ,t b \x ai t a ) = Xb VxVpe^ d < ipi - S ^). (2) 

J x (t a )=x a 

The new Hamiltonian H typically depends on first few momenta of u)(v,t). 
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Under normal circumstances, the smeared PI do not saturate CKE for 
Markovian processes, in fact, ad hoc choices of smearing distributions typ- 
ically introduce memory into a dynamics. Our goal here will be to isolate 
the general class of continuous smearing distributions that conserve CKE. 
In Ref. 7 we have shown that ensuing distributions have important applica- 
tions in evaluations of PI or in simplifications of the associated stochastic 
differential equations. In this note we shall briefly present the latter ap- 
plication. On the way we also mention a simple implication in economical 
models with stochastic volatility. 



2. Most general class of smearing distributions 

We look for uj(v,t) fulfilling CKE for any intermediate time t c : 

/oo 
dx P(x b , t b \x, t c )P(x, t c \x a , t a ) . (3) 
-oo 

It can be shown 7 that Eq.(3) is fulfilled only when the integral equation 

J dz 1 u(z',t)au(a(z - z'),^j = bw(bz,^j , (4) 
(a, b € R + and 1 + 1/a = 1/6) holds. By defining the Laplace image function 



uj as 



r-OO 

u)(£,t)= dz e - iz uj(z,t), 3?C>0, (5) 
Jo 

Eq.(4) can be equivalently formulated as the functional equation for Cj 

w&Q&feA =&(i t i) ■ (6) 



a' a J \ b ' b J 

Assumed normalizability and positivity of u)(v,t) implies that smearing 
distributions are always Laplace transformable. After setting a = l/a we 
get 

u>(£, t) Cj{ai, at) = w(£ + a£, t + at) . (7) 

This equation can be solved by iterations. An explicit general solution of 
Eq.(7) was found in Ref. 7 and it reads 

<~,(ft\ / [GKA)]*, when tjtO, 

- \««, when t = 0. (8) 
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G{x) is an arbitrary continuous function of x. Constant k is determined 
through the initial-time value of u)(v,t). In particular, t — > solution (cf., 
Eq.(8)) gives 

lim u(v,t) = 8(v + log k)5(v + log k) = S + (v + logK). (9) 

Let us also notice that (7) implies r) > for all r and £ which gives 
G(x) > for all x. This allows us to write 

[G(C/t)] t =e- F ^^ T . (10) 

Here F(x) is some continuous function of x. Final co(v, t) can be obtained via 
real Laplace's inverse transformation known as Post's inversion formula: 6 



, , ,. (-l) fe fk\ k+1 d k u>(x,t) 
w{v,t) = lim v ' 1 > 



k\ V v ) dx k 



(11) 

x—k/v 



Real inverse transform (11) is essential because the solution of the functional 
equation (7) was found only for real variables in Q(£,i). In fact, complex 
functional equations are notoriously difficult to solve. 

We finally point out that the result (8) is true also in the case when 
vH(p, x) — vHi(p, x) + H 2 {p 1 x), such that [H\,H2\ = 0. 

3. Explicit form of H 

To find H we use Post's formula (11) which directly gives 

P{Xb, tfr, X a: t a ^j 

■ Jo (fc-l)l ».,■>■ J x(ta)=x ' 



k^ooJo ' {k-l)\ dx k J x {t a )=x a 
fx(t b )=x b r roo 

= j VxVp j Ay u(y,t) e- y f*a drH 

Jx(t a )=x a YJo 



e l It ^ dT P* 



= / VxVp eti dT Vr*- F W) . ( i2 ) 

J x{t a )—x a 

In passing from second to third line we have used the asymptotic expansion 
of the modified Bessel function K k {\fkx) for large k. On the last line we 
have utilized the definition of the Laplace image. Result (12) thus allows to 
identify H with F(H). Let us finally mention that the normalization 

/>oo 

1=1 dv uj(v,t) =u)(0,i) (13) 
Jo 

implies that F(0) = 0. 
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4. Kramers-Moyal expansion for u>(v,t) 

Let us now mention a simple application in stochastic processes. To this 
end we notice that both P(xb,tb\x a ,t a ) and 

P v (x b ,t b \x a ,t a ) = f X(tb) X VxVp ett dr(i P x-„fO (u) 

J x(t a )—x a 

fulfil CKE, so they can be alternatively evaluated by solving Kramers- 
Moyal's (KM) equations: 1 - 8 

d 

—P(x b ,t b \x a ,t a ) =h K M P{xb,t b \x a ,t a ) , (15) 

Otb 

d 

-^—P v (Xbi tb\x a , t a ) = L^- M P v (xb, tb\x a , to) , (16) 
Otb 

where the KM operator L^m is 

L ^ = E(-^j D^{x b ,tb), (17) 

and similarly for ^ V KM - The coefficients and are defined through 
the corresponding short-time transitional probabilities, so e.g., D^\x,t) is 

1 1 f°° 

Di n \x,t) = -flim- / dy(y-x) n P v (y,t + r\x,t). (18) 
n! r^o r 7_ oc 

Equations (16) can be cast into an equivalent (and more convenient) system 
of equations. For this we rewrite (4) as 

u>{z,t)= &z' Lo(z\i!)P"(z,t\z' ,t'), (19) 
Jo 

with the conditional probability 
P»{z, t\z', t') = 6(tz - t'z') co (^f,t t> 



Jo ' 



dzP UJ {z,t\z',t') = l, lim P" (z, t + r\z',t) = 5+ (z-z'). (20) 



Eqs. (19)-(20) ensure that the transition probability P w (z, t\z', t') obeys 
CKE for a Markovian process. Since the process is Markovian, we can define 
KM coefficients in the usual way as 

1 r°° 

K^ n \v,t) = lim^y dx(x-v) n P w (x,t + T\v,t). (21) 
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From ensuing CKE for a short-time transition probability one may directly 
write down the KM equation for uj(v,t ba ) 

^ oo / r, \ n 

O , . s „ (UJ) , . s „ ful / O \ 

w(t 



dt 



ba 



(v,t ba )=^ M uj(v,t ba ), L^ = E(-|;) * (n) (»>*6a). (22) 

n— 1 ^ ' 



In cases when both (16) and (22) are naturally or artificially truncated at 
n = 2 one gets two coupled Fokkcr-Planck equations 



— u>(v,t) = jJ~plu}(v,t), —P v (x b ,t b \x a ,t a )='L FP P v (x b ,t b \xa,t a ), (23) 
pp = l^ K ^ M (n = 1,2), L^p = \^ KM {n = 1,2). 



On the level of sample paths the system (23) is represented by two coupled 
Ito's stochastic differential equations 

dx b = L>« {x b , t b ) dt b + \j 2D [ v ] {x b ,t b ) &Wi , 
dv = K^(v,t ba )dt ba + ^2K^){v,t ba )dW 2 . (24) 
Here W\{t b ) and W2(t ba ) are respective Wiener processes. 

5. Economical models with stochastic volatility 

Note that in (24) the dynamics of the variance v is explicitly separated 
from the dynamics of x b . This is a desirable starting point, for instance, in 
option pricing models. 1 As a simple illustration we discuss the stochastic 
volatility model presented in Ref. 2 To this end take G to be 

G{x)= i^Tl)' b€R+;ceR+. (25) 

This gives 

,<oo=(^y * »M)-fi^.-~. <-» 

F(0) = as it should. Distribution (26) corresponds to the Gamma distri- 
bution 2 ' 9 fu,ct{v)- The Hamiltonian H associated with (26) reads 

H{p,x)=vb\og(p^ + l}, (27) 

where v — c/b is the mean of u(v,t). As H we use the Hamiltonian from 
Refs. 1 ' 2 which has the form p 2 /2 + ip(r/w— 1/2), r is a constant. This choice 
ensures that P v (x b , t b \x a , t a ) represents a riskfree martingale distribution. 1 
Full discussion of this model without truncation is presented in Ref. 2 
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We consider here the truncated-level description that is epitomized by 
Ito's stochastic equations (24). The corresponding drift and diffusion coef- 
ficients D ( v 1} and Di 2) are then (cf. Eq.(18)) 

D£\x,t b ) = (r-|), D™(x,t b ) = \. (28) 

For the coefficients an explicit computation gives 

K^\v,tb a ) = j—{v — v), K^(v,t ba ) = -J--^. (29) 

Consequently, the Ito's system takes the form 

dx b =( r -%) dt b + VtJdWi, dv = ^-(v-v) dt ba + ^-\hr dW 2 . (30) 
v l) t ba t ba V o 

Let us now view x b as a logarithm of a stock price S, and v and r as the 
corresponding variance and drift. If, in addition, we replace for large t ab the 
quantity with \/v, the systems (30) reduces to 

dS = rSdt b +y/vSdWi, dv = 7 (v - v) dt ba + dW 2 . (31) 

The system of equations (30) corresponds to Heston's stochastic volatility 
model, 1,10 with mean-reversion speed 7 = l/t ba and volatility of volatility 
e = l/(t b aVb). 
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